Fluxon Pinning Through Interaction with the Superconducting Wiring of Long 

Annular Josephson Junctions 
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The statics and dynamics of magnetic flux quanta (fluxons) 
in long annular Josephson tunnel junctions have been inves- 
tigated. Pinning by interaction of the fluxon field outside the 
junction with the superconducting wiring has been observed 
in spatially resolving measurements using low-temperature 
scanning electron microscopy. We were able to influence the 
characteristics of this field by carefully modifying the beam- 
induced fluxon trapping procedure. In this way we were able 
to select the pinning site acting on the fluxon. 
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I. INTRODUCTION 

The nonlinear system constituted by a Josephson tun- 
nel junction has served as a rich source of interesting 
physical phenomena over the past decades jl| . Especially 
long junctions, i.e. junctions with at least one dimension 
large against the Josephson penetration depth Aj, have 
proven to be a suitable tool for the investigation of soliton 
(fluxon) statics and dynamics. Most of the experimen- 
tal work was carried out on linear Josephson junctions 
in the inline or overlap geometry The samples used 
for the experiments described in this article were of an- 
nular quasi-one-dimensional shape and are characterized 
in detail in section II. 

One-dimensional long Josephson junctions are gov- 
erned by the perturbed sine-Gordon equation ||^: 

y^xx - ftt = sin(^ + aifit - (iipxxt - 7 + f{^^ v), (1) 

where ip{x,t) is the spatially and time dependent phase 
difference across the tunnel barrier between the two 
Cooper pair systems of the superconducting electrodes. 
In the case of annular Josephson junctions, x runs az- 
imuthally around the junction. The variables x and t 
are normalized to the Josephson penetration depth Xj 
and the inverse of the plasma frequency uj~^, respec- 
tively. The dissipation parameters a and (3 account for 
losses due to quasi-particle currents across and parallel to 



the tunnel barrier, respectively, whereas 7 represents the 
normalized bias current. The function f{x, tp) describes 
inhomogenities of various origins, such as externally ap- 
plied magnetic fields or local changes of the thickness 
of the insulating barrier (e.g. microshorts). For an- 
nular junctions the dynamics and statics of solitons in 
the presence of inhomogenities have already been inves- 
tigated experimentally and numerically. Because 
strictly speaking the term 'soliton' only refers to certain 
exact solutions of the unperturbed sine-Gordon equation, 
we use the word 'fluxon' from this point on. In the an- 
nular geometry the boundary conditions for Eq.(|l|) are 
periodic, i.e. 



^(l,t) = (p{0,t) + n2Tr 



ip,{l,t) = ip.,iO,t), (2) 



where I is the ring circumference in units of Xj and n is 
the net number of fluxons trapped in the junction. This 
means that there can only be an integer number of flux- 
ons inside the junction at any given time. If each of the 
two superconducting junction electrodes is thicker than 
the magnetic penetration depth Xl (which is the case 
for the samples studied here), the total magnetic flux 
and thus the net fluxon number in the annular junction 
is conserved. This conservation and the fact, that flux- 
ons can move inside the junctions without influence of 
boundaries have motivated the experiments with annu- 
lar junctions. 

Even though there have been investigations to deter- 
mine the behavior of fluxons in externally applied mag- 
netic fields, there have been no experiments so far, to our 
knowledge, considering the fluxon magnetic field outside 
the actual junction. In our experiments we found ev- 
idence that an interaction of this field with the super- 
conducting wiring takes place and results in a pinning 
potential when junctions of the "Lyngby geometry" ||] 
are used, as is the case in most of the studies in the liter- 
ature. We have also investigated the single-fluxon statics 
in this potential by means of low temperature scanning 
electron microscopy (LTSEM) which allows a spatially 
resolved measurement of the Josephson current. These 
results are presented in section III. Moreover, we discov- 
ered consequences of this pinning in the dynamics of mul- 
tifluxon systems to the degree that a consecutive transi- 
tion of the fluxons into the dynamic state was observed 
as described in section IV. 
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II. THE SAMPLES AND THE EXPERIMENTAL 
SETUP 

The samples that have been used for our measurements 
are Nb/AlOx/Nb tunnel junctions with the annular 
geometry introduced by the Lyngby group [see Fig. Q [^ . 
The current density jc at 4.2 K of these junctions is 1000 
A/cm^, the mean radius is 90 /^m. The width of the rings 
is 10 /im (sample A) and 5 /im (sample B), respectively. 
The circumference of the ring, being the junction length, 
is 565 /im for both samples, corresponding to about 49 
Xj. This is very long compared to most other experi- 
ments conducted with annular junctions. This fact turns 
out to be vital for the observation of the phenomena de- 
scribed in this article. At 4.2 K the critical currents of 
the samples are approximately 45 mA (sample A) and 23 
mA (sample B). The zero field step asymptotic voltages 
were measured to be about 30 /iV at 4.2 K 

We performed spatially resolved measurements by 
means of low-temperature scanning electron microscopy 
(LTSEM) LTSEM allows the local thermal per- 

turbation of the junction due to the electron beam with 
the focus at xq during operation of the sample at liq- 
uid helium temperatures. The spatial extension of this 
perturbation is given by the thermal healing length and 
its value determines the spatial resolution, being about 
1-2 /im for the samples studied here. The beam-induced 
temperature increment AT(a;o) can be tuned by the elec- 
tron beam power. Since the thermal relaxation time for 
the beam induced local thermal perturbation is about 
100 ns, the sample response signal is a time averaged 
information about the junction dynamics. The latter 
evolve on a time scale of about 10 ps. 

In this article we present two different imaging tech- 
niques for the pinned fluxons. The first method uses the 
electron beam induced additional energy loss during the 
collision of two fluxons. In case of fluxon-antifluxon col- 
lisions a significant electron beam induced voltage signal 
AV{xq) < of the current-biased junction at the colli- 
sion sites can be observed [ |l^ , p^ . Notice, that we are 
dealing with the collision of fluxons of the same polarity 
in contrast to the fluxon-antifluxon collisions considered 
in Ref. [|l2|.^. Our experimental results show that also 
in the case of the collision of unipolar fluxons, the local 
thermal perturbation due to the electron beam causes a 
significant sample response AV{xo) < at the collision 
sites, as it has been already observed by Keil et al. [|j 

The second method is based on the two-dimensional 
imaging of the spatial distribution of the maximum dc 
Josephson current density p^ , p^ 

^ jc{x) sin <f{x) (3) 

depending on jc and on the phase difference ip. In our 
case, where jc does not show any signiflcant spatial de- 
pendence across the junction |16[| , j{x) is directly pro- 
portional to sin 1^9(2;). The local thermal perturbation 
induced by the electron beam at xq decreases jdxo) and, 



neglecting any nonlocal effect, the change Aj{xo) is a 
direct measure of smip{xo) |14[| . The junction is biased 
close to the total critical current 

Ic Ri I — / jcSinip{x) dx (4) 
Jo 

and Ic is continuously measured during scanning the e- 
beam across the sample. We obtain the electron-beam- 
induced signal 

- A/c(xo) cx smip{xo), (5) 

if the area perturbed by the beam is small compared to 
Xj . In case of a hysteretic current-voltage-characteristic 
{IV C) this imaging technique is described in [ p[ , 
whereas for a nonhysteretic IVC a description of the 
imaging technique can be found in [ p^ . For all images 
shown in this article any nonlocal effect of the local per- 
turbation due to the macroscopic quantum properties of 
a Josephson junction can be ruled out [0. For a com- 
prehensive review of the imaging of Josephson junction 
dynamics see Ref. [p^ . 

III. SINGLE-FLUXON EXPERIMENTS 

The key prerequisite for studying fluxon motion in an- 
nular Josephson junctions is fluxon trapping. For this 
purpose one has to break superconductivity. In addition, 
a magnetic field has to be applied during cooling of the 
sample through Tc- One or more fluxons can be trapped 
during this procedure. There are several methods for 
fluxon trapping p9|j20[ |. Since in a scanning electron mi- 
croscope wc investigate the cold samples, a reliable way 
of introducing fluxons into the junction consists of locally 
heating the superconducting electrodes to a temperature 
above Tc by electron beam irradiation. A magnetic field 
can be applied during cooling by passing some current 
through the junction or by means of an external solenoid. 
In our experiments we found it favorable to apply an ex- 
ternal field to achieve controllable trapping of single flux- 
ons. Figure |l| shows schematically the procedure of in- 
troducing a fluxon into the junction. After the scanning 
unit of the LTSEM is adjusted to deflect the electrons 
onto one of the electrodes the beam is switched on. This 
operation guarantees that the area heated above Tc by 
the beam does not contain any magnetic flux because it 
is generated well apart from the edges in the supercon- 
ducting material. Subsequently, when the electron beam 
is moved towards the center of the ring, the second elec- 
trode is heated as well. In this electrode magnetic flux 
supplied by the external field can be dragged along. Note 
that the different starting points depicted in Fig. |l] (a) 
and (b) lead to different configurations of the external 
fiuxon field. Inside the junction the results of the two 
procedures are identical. 

After fluxon trapping, the number of trapped fluxons 
can easily be determined experimentally by recording the 
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IVC of the junction through applying a current to the 
junction. Due to the Lorentz force between the bias cur- 
rent and the magnetic moment of the fluxon the latter is 
accelerated, moves through the junction, and a voltage 
drop is observed. The IVC of sample A with one fluxon 
trapped by the procedure described above is shown in 
Fig. ^ (a) as well as the differential resistance (dU /dl) 
measured through electronic differentiation with a lock- 
in amplifier. The asymptotic voltage of about 30 fiV 
corresponds to the limiting Swihart velocity c |2^, the 
velocity of light in the junction transmission line. Note 
that the curvature of the IVC does not comply with 
perturbation theory, which predicts d^U/dP < f for all 
V. We believe that the reason for this is the voltage de- 
pendency of the quasiparticle tunneling probability [^2| , 
which means that the dissipation parameter a is voltage 
dependent as well. 

The small critical current of about 1 mA (roughly 3 
% of Ic without fluxons trapped) suggests the existence 
of a pinning potential. Various authors have reported 
similar observations [p|p^-[25|]. Davidson et al. found a 
small critical current for the IVC of a junction with a 
length of 15 Xj when fluxons were trapped. That critical 
current depended on the junction history 

The statics of a single fluxon in the pinning potential 
responsible for the small critical current described above 
was investigated by LTSEM. In Fig. | (b) and (c) LT- 
SEM images of a single fluxon trapped by the procedure 
depicted in Fig. |l] (a) are shown (sample A). They show 
a single fluxon at two different positions in the ring. The 
white areas correspond to a Josephson supercurrent into 
the paper plane, the black areas into the opposite direc- 
tion. The magnetic moment and the direction in which 
the Lorentz force acts are symbolized by arrows. Appar- 
ently the location of the pinning center depends on the 
direction of the bias current. Considering that, because 
of the way the fluxon was trapped, the magnetic field of 
the fluxon encloses only the upper electrode (see Fig. |l| 
(a)). Therefore, we conclude that the fluxon gets pinned 
through the interaction with the upper electrode. 

Further evidence for the validity of this argument is 
given by the LTSEM images of Fig. ^ The fluxon shown 
in these images was trapped in the junction using the 
procedure depicted in Fig. |l| (b) (sample A). Hence the 
magnetic field of the fluxon encloses only the lower elec- 
trode and gets caught on that wiring when acted on by 
the bias current. The fact, that the position of the pinned 
fluxons with respect to the horizontal symmetry axis of 
the junction differs for reverse bias current directions in 
both Figures ^ and ^, hints at the existence of additional 
forces on the fluxon. The measurements took place in 
a reasonably well shielded environment. However it is 
possible, that some external magnetic field is present at 
the sample site, which would act as an additional pinning 
force. 

An investigation of single fluxons trapped in the junc- 
tion by crossing both electrodes from the side, i.e. switch- 
ing on the electron beam outside the electrodes and then 



moving it into the center of the ring by crossing the 
narrowest part of the combined electrodes, confirms the 
above interpretation. This procedure allows the field 
of the fluxon to circle both electrodes as schematically 
shown in Fig. ^. Even though shielding currents flow in 
the electrodes to ensure the unambiguity of the Cooper 
pair function in each electrode (the magnetic flux con- 
tained by a single fluxon is exactly the elementary flux 
quantum <f>o), pinning is still expected to occur. The rea- 
son is the fluxon field outside the junction in the shape 
of a distorted magnetic dipole field having its maximum 
value just at the outer border of the junction (see Fig. 
I (b)). The distortion is due to the Meissner effect in 
the superconducting electrodes. Inside the junction the 
fluxon field is localized to an azimuthal length compa- 
rable to the Josephson penetration depth. In this case 
interaction of the fluxon magnetic field with the wiring 
of both electrodes takes place. The resulting pinning 
potential is expected to be shallower than the potential 
discussed earlier for the procedures in Fig. ^ (a) and (b). 
Consequently the critical current should be lower. Ex- 
periments showed indeed, that trapping of single fluxons 
over the sides resulted in a value of Ic being 50 % smaller 
than for trapping from the middle of the electrodes. 

It is interesting that the transition into the dynamic 
state does not happen in an abrupt manner. Even if a 
significant voltage drop is measured across the junction 
for a single-fluxon mode (for example a voltage greater 
than 5 fiV in Fig. | (a)) the LTSEM images (not shown 
here) still display the Josephson current distribution of a 
fluxon at rest. We believe that this is a result of the rel- 
atively weak pinning forces which allow the fluxon to get 
depinned only part of the time at first. As described in 
section II the LTSEM images are time averaged measure- 
ments. Therefore, the images recorded at finite voltages 
can be interpreted as time averaged images of a sporad- 
ically moving fluxon. 



IV. MULTIFLUXON EXPERIMENTS 

The features of IV Cs do not change qualitatively when 
more than one fluxon is trapped. Figure ^ (b) shows 
the IVC and the differential resistance for the case of 
two fluxons trapped by the procedure shown in Fig || 
(a) . The dU/ dl-cmve displays two peaks at low voltages 
corresponding to the consecutive transition of the two 
fluxons from the static into the dynamic state. Note that 
Ic of the two fluxon system is only about 60 % of the 
value for one fluxon. It appears that this is due to the 
interfluxon repulsion of the fluxon magnetic moments. 

i- 

In the light of the concentration of the extrajunctional 
fluxon field in the immediate vicinity of the junction (see 
Fig. H (b)) it is noteworthy that it was not possible to 
trap more than seven fluxons in the junction by repeat- 
edly applying the procedures shown in Fig |^. On the 
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other hand it needed only little effort to trap up to 50 
fluxons by introducing them over the sides i.e. allowing 
the fluxon magnetic field to enclose both electrodes . 
As already explained in section 111 (reduced) pinning is 
still expected for this case. The IVC and the dU/dl- 
curve for such a multifluxon mode show a rather compli- 
cated shape as depicted in Fig. ^ for the situation of sev- 
enteen fluxons trapped in the junction (sample B). Due 
to the presence of a small ohmic resistance in series with 
the junction in this particular measurement the IVC is 
tilted and the asymptotic voltage is larger than the 510 
/iV expected for a seventeen fluxon system. Nevertheless 
the structure of the dU /d/-curve can be explained by the 
consecutive transition of all seventeen fluxons from the 
static into the dynamic state. 

Fluxon-fluxon collisions are expected for the situation 
where more than one fluxon is trapped in the junction, 
while at least one of them is moving and the others are 
pinned. A large number of LTSEM images were recorded 
near the resonant structures in Fig. |6[ Arrows indicate 
the voltages at which selected images, shown in Fig. |^, 
were taken. In the LTSEM images we observe the follow- 
ing process. One after another the fluxons get depinned 
when the bias current is increased, as indicated by the 
decreasing number of signal peaks. Finally only a single 
one is still at rest while sixteen are moving as seen in Fig. 
(d) . The moving fluxons are used as a detector for the 
pinned ones by the collisions that take place once during 
each revolution of every moving fluxon. 

At this point it is necessary to recall that the tran- 
sition into the dynamic state is not abrupt. Instead a 
sporadic change between statics and dynamics of a par- 
ticular fluxon takes place. Therefore it makes sense that 
the image in Fig.|^ (d) still shows the last fluxon, even 
though it was recorded at a voltage greater than the one 
corresponding, to the last peak in the differential resis- 
tance in Fig. |g. This is in accordance with the findings 
for single trapped fluxons. It is interesting that the po- 
sition of the last static fluxon shown in Fig. (d) is 
identical with the one of a single fluxon getting pinned 
at the wiring of the upper electrode. For reversed bias 
current direction the pinning of the last fluxon at the 
other side of the same wiring was observed. 

Furthermore we want to point out the way in which the 
size of the fluxons is different in each of the images of Fig. 

. The fluxon compression (see eq. |^ results in a reduced 
fluxon length, since the fluxons are 'squeezed' into the 
junction. The more fluxons change to the dynamic state, 
the less compression the static ones experience. There- 
fore the static fluxons expand. Finally the last static 
fluxon (Fig. 1^ (d)) reaches the 'relaxed' size, i.e. the size 
we observe if only a single fluxon exists in the junction 
at rest, being approximately 27rAj. 



V. CONCLUSIONS 

We have investigated the statics and dynamics of 
(Josephson) fluxons in long annular Nb/AlOx/Nb tun- 
nel junctions. The current-voltage characteristics showed 
evidence of fluxon pinning. Applying low temperature 
scanning electron microscopy, we were able to image the 
statics of single pinned fluxons. By carefully control- 
ling the procedure through which fluxons were introduced 
into the junction, it was possible to influence the posi- 
tion of the pinned fluxons. These experiments show that 
this pinning is a result of the interaction of the mag- 
netic fluxon field with the superconducting wiring of the 
junction. We have also investigated the behavior of mul- 
tifiuxon states. For the case of 17 fiuxons trapped in the 
junction we found that the fluxons get depinned one after 
another with increasing bias current. 
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FIG. 3. (a) Sketch of the annular junction with the same 
orientation as in (b) and (c) but drawn to different scale. 
The voltage images in (b) and (c) show a single static fluxon 
trapped by the procedure of Fig 1 (a) (sample A). For imag- 
ing, the junction was biased at I = -1-1.9 mA and I = -2.2 mA, 
respectively. Positive bias current values denote a current flow 
from the top to the bottom electrode. Fluxon magnetic mo- 
ment and the direction of the Lorentz force are symbolized by 
arrows. 

FIG. 4. LTSEM voltage images of a single static fluxon 
trapped by the procedure depicted in Fig. 1 (b). The junc- 
tion (sample A) was biased at I = -1-1.7 mA (a) and I = -1.6 
mA (b), respectively. In contrast to the single-fluxon config- 
uration shown in Fig. 3 pinning at the edges of the lower 
electrodes can be observed (sample orientation as shown in 
Fig. 3 (a)). 

FIG. 5. (a) Schematic drawing of a possible extrajunctional 
field configuration for the case of a fluxon trapped by cross- 
ing both electrodes simultaneously with the focused electron 
beam at their narrowest point, i.e. along the vertical sym- 
metry axis in Fig. 3 (a). Shielding currents in the electrodes 
are necessary to ensure the unambiguity of the Cooper pair 
wave functions, which results in flux quantization in the rings 
of the electrodes, (b) Sketch of the flux distribution for the 
upper electrode in the situation depicted in (a). Even though 
the net flux in the opening $„et = 0, the fluxon magnetic fleld 
outside the junction will still interact with both wirings. 



FIG. 1. Schematic representation of the trapping proce- 
dures for fluxons resulting in a defined extrajunctional mag- 
netic field configuration. The electron beam is switched on 
when already focused on one of the electrodes. The energy 
deposited by the electron beam results in local heating to 
temperatures above Tc. Magnetic flux can only penetrate the 
normal region in the electrode crossed second by the beam. 

(a) shows the procedure starting from the bottom electrode, 

(b) from the top electrode. 

FIG. 2. Current-voltage curve {IVC) and differential re- 
sistance (dU/dl) of the annular junction (sample A) with (a) 
one fluxon and (b) two fluxons trapped using the trapping 
procedure of Fig. 1 (a). Small critical currents of about 1 
mA and 0.6 mA, respectively, (approx. 3 % of 7c with no 
fluxon) and the peaks in the dU/dl curves are evidence of 
fluxon pinning. T fs 4.2 K. 



FIG. 6. Current-voltage characteristic (IVC) and differen- 
tial resistance (dU/dl) of the junction (sample B) with 17 flux- 
ons trapped. Due to a small ohmic resistance in series with 
the junction thelVC is tilted. The peaks in the dU/dl-curve 
are a result of the consecutive transition of the fluxons from 
the static to the dynamic state. The arrows indicate the volt- 
ages at which the LTSEM images in Fig. 7 were recorded. 

FIG. 7. LTSEM voltage images of 17 fluxons trapped in the 

junction (sample B) at different bias current values, (sample 
orientation as shown in Fig. 3 (a)) (a) All 17 fluxons are 
at rest. The size of the fluxons is reduced compared to the 
single-fluxon case of Fig. 3 and 4 due to fluxon compression 
inside the junction, (b) At I = 2.1 mA the image still shows 
eleven signal peaks, indicating that six fluxons are moving, 
(c) Four peaks are left at a bias current I = 3.6 mA. (d). 
Finally at I = 7.3 mA only a single static fluxon is visible. 
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